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Abstract 

00 ' 

Q\ \ Using microscopic models in which both photons and excitons are treated as microscopic 

degrees of freedom, we discuss polaritons of two cases: One is the case when excitonic 
parameters are time dependent. The time dependence causes creation of polaritons from a 
Ph. "false vacuum." It is shown that both the creation sepctra and the creation efficiency are 

much different from the results of the previous studies. The other is polaritons in absorptive 

7— i \ and inhomogeneous cavities. A polariton in such a system cannot be viewed as a back and 

forth oscillation between a photon state and an exciton state. 
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1. Introduction 

o . 

Optical phenomena in condensed matter are usually discussed in the case when the optical 
field can be treated as a classical field. Recently, however, quantum optical phenomena 
in condensed matter, in which the optical field exhibits its quantum natures in condensed 
matter, have been attracting much attention. A simple theoretical framework, which is 
widely used in the literature, for the analysis of such phenomena is to somehow quantize the 
macroscopic electromagnetic (EM) fields in the "macroscopic Maxwell equations" (Maxwell 
equations in matter) Jl| . The simple theory treats the matter as an effective medium, whose 
. properties are assumed to be fully described by phenomenological parameters such as a 

dielectric constant and nonlinear susceptibilities. The theory works well for simple problems. 
However, its validity is not clear when dispersion and absorption is important, when the 
system is inhomogeneous, when nonlinearities are important, or when the phenomenological 
parameters are time dependent. 

Another way of formulating quantum optics in condensed matter is a microscopic approach. 
In this case the matter is represented as microscopic polarization fields, and both the EM 
fields and the polarization fields are quantized. Upon diagonalization of the Hamiltonian of 
the coupled photon-polarization system, Hopfield obtained "polaritons" as eigenstates of 
the total system. His theory has been extended by many researchers to study polaritons of 
various types. For example, polaritons in the case when excitons are subject to dissipation 
were studied by Huttner, Baumberg and Barnett ||. 

In this work, we discuss two cases: One is polaritons in the case when excitonic parameters 
(the exciton energy and exciton-photon coupling) are time dependent |1J. The time depen- 
dence causes creation of polaritons from a "false vacuum." This phenomenon is analogous to 
the dynamical Casimir effect, which is creation of photons by quick movement of mirrors in 
vacuum 0, || . We will present the creation spectra of both the lower- and upper-branch po- 
laritons. The other is polaritons in absorptive and inhomogeneous cavities . We point out 
that a polariton in such a system cannot be viewed as a back and forth oscillation between 
a photon state and an exciton state, and a photon state will evolve into a superposition of 
many, different exciton states, and will not revive. 



2. Creation of polaritons from a false vacuum of solids 



It is usually assumed in discussions of polaritons that parameters such as the exciton 
energy or exciton-photon coupling are constant, i.e., independent of time. What happens 
if the parameters are time dependent? We can show that polaritons will be created even 
when the initial state is the vacuum. This is a general property of quantum theory: Particles 
will be created when the Lagrangian (or, almost equivalently, equations of motion) has a 
time-dependent parameter(s) f|. This can be understood simply as follows. Suppose that 
the parameters are constant for t < 0. For this constant values of parameters we can (in 
principle) find out the ground state, which is called a "vacuum," in which no polaritons 
are present. This vacuum is denoted by |0i). Assume that the parameters vary during 
< t < T, and again become constant for t > T. Then for t > T we can find out another 
ground state, which is denoted by |0 2 ). Under the initial condition that the state vector \ip) 
is 1 1 ) for t < 0, we may solve the Schrodinger equation ||. Then we will find that \ip) at 
t > T differs from |0 2 ) in general. This means that we have a finite probability of finding 
particles at t > T. More particles are normally created for faster variations of parameters 
(i.e., for shorter T). This is understood if we think of the limiting case of T — > oo, for which 
|^) evolves adiabatically from |0i) at t < to [O2) at t > T, and no particles are created. 
That is, non-adiabatic (fast) variation is necessary for particle creation. 

We here consider the particle creation in condensed matter. Such a phenomenon is known 
in quantum optics as "spontaneous parametric fluorescence" 0. The conventional theory of 
the spontaneous parametric fluorescence is a phenomenological one. That is, it relies on the 
phenomenological quantization scheme, assuming that dispersion and absorption are absent. 
It also assumes a phenomenological interaction Hamiltonian of the following form; 

H int = galaldp + h.c, (1) 

where a s , a, and a p are annihilation operators of "signal" (s), "idler" (i) and pump photons, 
respectively, and g is an effective coupling constant which is proportional to the second-order 
nonlinear susceptibility. Through this interaction, the equations of motion of the signal- and 
idler-photon fields are modulated at the frequency u p (= u s + Ui) of the pump photon 
field. That is, the time scale T of the modulation is T ~ 1/uj s ~ 1/W This is a very 
fast modulation, and the signal and idler photons can be created efficiently enough to be 
observed experimentally. 
Similar situations have been studied by Yablonivitch [ID and by Schwinger [IT|. They 



investigated photon creation in the cases when the dielectric constant e of material is time 
dependent. The time dependence may be due to creation of an electron- hole plasma by a 



short laser pulse [ 1C ] , or by sudden collapse of a bubble in a liquid [11]. As in the case of the 
spontaneous parametric fluorescence, the theoretical analyses are phenomenological ones. 
That is, they rely on the phenomenological quantization scheme, assuming that dispersion 
and absorption are absent. It was furthermore assumed that e changes either discontinuously 
(T = 0) [1TDL [TT], [T2|, or almost adiabatically (T ~ 00) |TUJ. The discontinuous change is 
unphysical, and, as will be shown below, turns out to lead to enormous overestimate of the 
particle-creation efficiency 



Because of these many assumptions, the previous theories ||, [TI| [Il|] have only a limited 
range of applicability. For example, they cannot be applied to the case when relevant photon 
energies hu are close to the excitonic energy hui x , because the dispersion is strong for uj ~ uj x . 
However, such a region of the photon energy should be most interesting. Such an interesting 



which is described by the Lagrangian density of the system of the exciton field X(r, £), the 
photon field A(r, t), and the longitudinal field U(r,t): 
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where the electric and magnetic fields are given by 

E = —A — VU, B = V x A. 



(2) 



(3) 



This Lagrangian is a generalization to the time-dependent uj x and a of the model proposed 
by Huttner et al. JJ. Following them, we apply the standard quantization procedure as 
follows. We impose the Coulomb gauge condition, so that A becomes transversal and U can 
be eliminated. By decomposing X into the transversal and longitudinal parts, 

X = X T + X 11 , (4) 

we can decompose the Lagrangian L into the transversal and longitudinal parts: 

L - I Cdv = L T + L" . (5) 



Since A is involved only in L T , we focus on this part. In the Coulomb gauge, it is convenient 
to quantize the fields in the k space: 



A(r, t)=Y, A(k, t)e ikr , where A(k, t) = A f (-k, t), 

k 

and similarly for X T . The fields are further decomposed into two polarizations; 

A(k, t) = £ A x (k, f)e A (k), k • e A (k) = 0, 

A=l,2 



(6) 



(7) 



and similarly for X T . By expressing L T in terms of these components, we can find their 
conjugate momenta U x (k,t) and P x (k,t) by differentiating L T with respect to A x (k,t) and 
X x (h,t), respectively. The Hamiltonian H (for the transversal parts) is obtained as 



H = ^(n A i A + P X X X ) - L T . 



(8) 



kA 



Note that this H, like L, has explicit t dependence, H = H(t), through the t dependent 
parameters. The fields are then quantized, in the standard manner, by imposing the equal- 
time commutation relations; 



[A x (k,t),U x \k',t)] = [X A (k,t),P A '(k',i)] =ihS x , x ,S kiV 
This leads to the Heisenberg equations of motion; 

ih^A x = [A x , H(t)], zh^X x = [X\ H(t)]. 



(9) 



(10) 



So far, the calculations are parallel to those of Ref. 0. The peculiar features of time- 



each t. Let a A (k, t) and & A (k, t) be the operators which diagonalize the Hamiltonian at time 
t: 

H{t) = hJ2 [ti a (k,t)a M (k,t)a x (k,t) + Cu b (k,t)b M (k,t)b x (k,t)} + c-numbers. (11) 

k,A 

If a and u x were constant, a x and b x would have become the usual annihilation operators 
(times a phase factor, which of course is unimportant) of lower-branch (LB) and upper- 
branch (UB) polaritons, respectively, with u a and uuj, being their eigen-frequencies. However, 
this is not the case because our a and u x depend on t. Moreover, it should also be emphasized 
that the equations, 

ih^a x = [5 A ,#00], thp x = [b x ,H(t)], (wrong) (12) 

which look like Eq. (|i~0[), does not hold for a x and b x . For these reasons, the diagonalization 
of H does not solve the problem. 

To solve the problem, we must first specify the initial state vector. For this purpose, we 
assume that the parameters are constant in the remote past: 

a(t) ~ oti, u) x {t) ~ uJi as t — > — oo. (13) 

This ensures that as t — > — oo the operators a x and b x and the frequencies u a and Ub approach 
the annihilation operators a x and b x and the eigen frequencies £o al and u bl , respectively, of 
the usual polariton operators, which are for the values of the parameters a — a.\ and uj x = uj\. 
We consider the case when the initial state vector is the vacuum of these polaritons: 

\ip) = |0i), where S A |0i) = 6 A |0i) = 0. (14) 

The state vector does not evolve because we are working in the Heisenberg picture. On the 
other hand, to specify the final states (the vacuum, one-particle states, two-particle states, 
etc.), we further assume that the parameters are constant also in the remote future: 

a(t) ~ «2, &x(t) ~ 0J2 as t — > +00. (15) 

Hence, as t — > +00 the operators a x and b x and the frequencies Cj a and Ub approach the 
annihilation operators a x and 6 A and the eigen frequencies u a2 and tD fe2 , respectively, of 
another polariton operators, which are for the values of the parameters a = a 2 and u x = uo 2 . 

We are interested in the number of created polaritons of each polariton mode in the final 
state. They are given by 

n A 2 (k) = (Oil^CkJaJCkJIOx), n x b2 (k) = (Oil&^k^k)^), (16) 

for each mode of the LB and UB polaritons, respectively. To evaluate these numbers, we 
must express a 2 and b 2 in terms of a x , b\ and their Hermitian conjugates. For our time 
dependent Lagrangian this relation takes the form of a Bogoliubov transformation, hence 
n a2 and hyi become finite 0. The coefficients of the transformation can be found by solving 
the Heisenberg equations of motion for the fields, Eq. (|10D, and inserting the results into the 
relation between the fields and a A , 6 A , a x ^ and b x ^ . 



We have performed a numerical calculation taking the following forms for u x (t) and a(t); 



(f <0) 

[^-sin(^)] (0<t<T) 



rl j " " * w 2 (t>T) { U> 

«i (t < 0) 

a{t) = { ai -M[M_ sin (^)] ( <t<T). (18) 
a 2 (t > T) 

These are continuous up to the second derivatives, whereas the third derivatives are discon- 
tinuous at t — and T. Figure 1(a) shows the calculated results for the number of created 
polaritons of each mode of the lower branch, n^ 2 (k), for various values of U\T. Figure 1(b) 
shows those of the upper branch, fi^ 2 (k). In both cases, strong excitonic features appear at 
ck/uji ~ 1 (which also means ck/u^ ~ 1 because we have taken uj\ ~ u 2 ). When ck < uoi 
(ck > uji), UB (LB) polaritons are created more efficiently than LB (UB) polaritons. This 
crossover occurs because we are here varying the exciton parameters, and the UB polariton 
has more (less) exciton component for ck < Ui {ck > Uj). These characteristics are correctly 
described only by a microscopic model. 

We also find that n^ 2 decreases very quickly as T is increased. To understand the decrease, 
let us investigate the case when k = 20ui/c and ai = a 2 = O.OOlcui^/eop. For such large 
k and small a, a LB polariton is almost an exciton, and we can understand the physics 
clearly. The dotted line of Fig. 2 plots n* 2 for this function of lo{F . We find that 

n^ 2 decreases exponentially as T is increased. We also find small oscillations for large T. 
This is found to be due to the fact that our u x (t) has singularities (third-order derivative is 
discontinuous at t — and T). To verify this we also consider the case when u x (t) takes the 
following analytic form 

<•*<«) - <•* - T^Wy < 19 > 

where T = lOr. (The coefficient is chosen in such a way that a small-T behavior agrees with 
the case of the non-analytic uj x .) To simplify discussion, let us take a = 0, for which a LB 
polariton becomes a pure exciton, and n^ 2 becomes independent of k. We then obtain the 
exact solution: 

- x = ^ 2 KK-"2)r] 

sinh|27rci;iT| sinh|27rc<j2Tj 

This is plotted in Fig. 2 with solid line. It is seen that the non-analyticity of the form 
(|17D induces oscillations. It was also shown |§ that for large oj\T (> 20) the non-analyticity 
enhances n\ 2 by many orders of magnitude. Since non-analytic change of u x (t) is unphysical, 
the result of analytic change of uj x (t) should be more convincing. However, most previous 
studies [|1(], [TTj, [L2| assumed that a time-dependent parameter (such as the dielectric constant) 
itself has discontinuities . Our results strongly suggest that such strong singularities would 
lead to physically incorrect results. Regarding our results of Fig. 1, they may be convincing 
because our slight singularity does not play important roles for such short T, as seen from 
Fig. 2. 

To summarize this section, we have considered a polariton system in which excitonic pa- 
rameters are time-dependent. The time dependence causes creation of polaritons from a 
"false vacuum," and we have evaluated the number of created polaritons for each polari- 



important, and the dielectric constant exhibits strong dispersions, we have employed a mi- 
croscopic model in which the polarization degrees of freedom are included as microscopic 
variables. Moreover, whereas most previous studies assumed sudden changes (T = 0) of a 
parameter (s), we have assumed more realistic situations in which excitonic parameters vary 
within a finite time (T > 0). Our results strikingly differ from the previous results, both 
qualitatively (strong excitonic features at ckju\ ~ 1) and quantitatively (by many orders of 
magnitude). 

3. Polaritons in absorptive and inhomogeneous cavities 

It is usually assumed in discussions of polaritons that the excitation level, which couples 
to photons, of the material is a discrete level. The excitation level may be an exciton level in 
the case of an excitonic polariton, or a phonon level for a phonon polariton. As the photon 
energy is increased to a continuous absorption spectrum, the single-level approximation 
breaks down. Even in such a case, however, we may treat, to a first approximation, the 
material excitations as bosons (polarization fields) because the excitations are composed 
of bosons or pair excitations of fermions. We may therefore write the Lagrangian density, 
assuming a one- dimensional system for simplicity, as 

£ = <CeM + £-mat + £int, (21) 

Cem = ?U 2 -c 2 ^Vl, (22) 




£mat = ~J q dup^X^^Xl), (23) 

C int = - / cLuauAXu, (24) 
Jo 

where A(x, t) is the vector potential of the electromagnetic (EM) field, and X u (x, t) denotes 
the polarization field of frequency u. Since we are considering the case when the photon 
energy lies in a continuous absorption spectrum, we have taken X u (x, 1) to have a continuous 
label u), and integral over u is performed in Eqs. (^) and (p4|). Moreover, we are interested 
in the case when the spatial distribution of the material is inhomogeneous because in such a 
case the material constitutes a lossy cavity, in which behavior of polaritons should be very 
interesting. When the material exists only in the regions of \x\ > £/2, we may express this 
inhomogeneous distribution by imposing 

X UJ (x,t) = for\x\<£/2. (25) 

The equal-time canonical quantization of A and X^ can be performed in the standard man- 
ner, and we have found that the Hamiltonian is diagonalized as [0 

H=Y,Hl duohuoa {qa \uo) ] a {qa \uo). (26) 

a=± q J 

Here, (u) is the annihilation operator of a polariton, which is given by 

j 



where a = ±1 and j = 0, 1, 2, • • • (or, a and q = 1, 2, • • •) label some modes of the cavity, 
—6oE and P are conjugate momenta of A and X, respectively, and (5, /3,7 and 7 are some 
coefficients. (For details of the notations, see Ref. 0.) It is seen that a^(u) is a superposi- 
tion of infinite number of modes (of A and X). In particular, the integration over continuous 
exciton modes results in the breakdown of the standard picture, that a polariton is a back 
and forth oscillation of photon — > exciton — > photon — > exciton —>■■■■. That is, when we 
prepare a photon state as an initial state, the state will evolve into a superposition of infinite 
number of exciton modes, and never returns to the initial photon state. 

Using the solution of such unusual polaritons, we can investigate various physical phenom- 
ena in absorptive and inhomogeneous cavities. For example, we have evaluated the radiative 
lifetime of an excited atom in such a cavity Calculations of other quantities are in 
progress. 
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Figure captions 

Fig. 1. The number of created polaritons per mode for (a) lower- and (b) upper-branch polaritons. 
The polariton parameters are taken as W2 = 0.95cji, a.% = 0.05ui^/eop, «2 = 0.03^1-^0/5, and 
ujiT = 0.1,2,4. 

Fig. 2. The number of created polaritons per mode for lower-branch (LB) polaritons is plotted 
as a function of u)\T. Dotted line is the result for the LB polariton of k = 20u;i/c, for a\ = 02 = 
O.OOlui-y/eop. Thin line represents the result when uj x (t) takes the analytic form, Eq. (|l^), and 
a = (which means that a LB polariton in this case is a pure exciton). In both cases 0J2 is taken 
to be 0.9wi. 
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